Evidence for a new particle on the worldsheet of the QCD flux tube 
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We propose a new approach for the calculation of the spectrum of excitations of QCD flux tubes. 
It relies on the fact that the worldsheet theory is integrable at low energies. With this approach, 
energy levels can be calculated for much shorter flux tubes than was previously possible, allowing 
for a quantitative comparison with existing lattice data. The improved theoretical control makes it 
manifest that existing lattice data provides strong evidence for a new pseudoscalar particle localized 
on the QCD fluxtube - the worldsheet axion. 



Recent advances in lattice simulations of quantum 
chromodynamics (QCD) have allowed to visualize confin- 
ing strings quite vividly 1 1 and to measure the spectrum 
of their low lying excitations with impressive accuracy [2] . 

These lattice results have lead to an embarrassing 
situation for theorists. On the one hand, even for 
strings whose length is merely twice their width, many 
of the energy levels in the lattice simulations show re- 
markable agreement with the energy levels of a the- 
ory that the QCD string is certainly not described by, 
the bosonic string defined at the quantum level through 
light-cone quantization [SJ. This Goddard-Goldstone- 
Rebbi-Thorn (GGRT) string is only Lorentz invariant in 
D = 26. QCD flux tubes, however, originate from a rela- 
tivistic theory in four dimensions and must be described 
by a worldsheet theory that respects Lorentz symmetry 
in D = 4. The agreement is thus rather surprising. On 
the other hand, existing theoretical techniques for cal- 
culating the flux tube spectra for the Lorentz-invariant 
Nambu-Goto (NG) string [IHB] break down for the rela- 
tively short strings that can be simulated with current 
lattice techniques. To make matters more confusing, 
there is also a family of energy levels which disagree badly 
with the predictions made by the GGRT theory [2]. 

This Letter describes a theoretical framework to com- 
pute energy levels of the NG string for much shorter 
lengths than previously possible. Our better theoretical 
understanding allows us to explain both why there was 
agreement between QCD flux tubes and GGRT strings 
for many levels, and why there was disagreement for oth- 
ers. Rather interestingly, we show that the data implies 
the existence of a massive pseudoscalar resonance on the 
string worldsheet. The NG string itself is thus insufficient 
to describe QCD strings. We explain how to include this 
resonance into our framework and measure its mass and 
width from the data. 

Before presenting our results, let us first describe the 
lattice data and compare it to the standard perturbative 
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FIG. 1: This plot shows AE = E- R/f s as a function of the 
length of the flux tube R for three different levels. The ground 
state is shown in green. States with spin 1 and with one and 
two units of longitudinal momentum are shown in orange and 
red, respectively. The value of £ 3 was determined from the 
ground state data. The solid line shows the prediction of a 
derivative expansion. The dashed lines shows the prediction 
of the GGRT theory. For the spin 1 states the prediction for 
a free theory is shown as dotted lines. 



results [1H5]. All numerical results we discuss are taken 
from [2] and are for gauge group SU(3). To measure 
the flux tube spectrum on the lattice, one calculates the 
discretized Yang-Mills 2 partition function with a Wil- 
son loop inserted at time r = and its conjugate at 
t = T, both winding around the compact spatial dimen- 
sion with periodicity R. From the asymptotics of the 



2 To avoid the string breaking due to quark pair production lattice 
1 See http://www.physics.adelaide.edu.au/theory/staff simulations are performed in the pure glue theory. We refer to 

/leinweber/VisualQCD/Nobel/ for animations. this theory as QCD. 
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partition function in the limit of large T, one deduces 
the energy of the ground state of a closed flux tube of 
length R. To measure the energies of excited states one 
deforms the shape of the Wilson line to project out lower- 
lying energy levels. The ground state energy as a func- 
tion of the string length R is shown in green in Fig. [T] 
The data agrees remarkably well with the GGRT ground 
state energy. Though surprising at first sight, this agree- 
ment finds a straightforward explanation in an effective 
field theory approach. The excitations of a long QCD 
flux tube are Goldstone particles. They arise because 
the presence of a long, straight string spontaneously 
breaks the target space Poincare symmetry 750(1, 3) to 
ISO(l, 1) x 50(2). For a recent discussion emphasizing 
this viewpoint, we refer the reader to [TJ. The standard 
method of calculating the effective string spectrum then 
is a derivative expansion, or cquivalently an expansion in 
the small parameter £ S /R. The non-linear ly realized tar- 
get space Lorentz symmetry imposes strong restrictions 
on the coefficients in this expansion and predicts all the 
coefficients up to £*/R 5 [5J. These universal coefficients 
in the expansion are the same as those in the expansion of 
the GGRT ground state energy. Figure [l] illustrates that 
this universality alone is enough to explain the ground 
state data. 

A bigger puzzle arises for excited states. Figure [T] 
shows the energy as a function of length for excited states 
with a single left-moving phonon with one and two units 
of longitudinal momentum. As before all the terms up 
to order £^/R 5 in the £ s /i?-expansion of the energies of 
these states are universal and coincide with the corre- 
sponding terms for the GGRT string. They agree well 
with the data at large R. However, for short strings the 
£ s /i?-expansion breaks down, and the universal terms no 
longer provide a good description of the data. 

By itself the breakdown of the perturbative expansion 
for small R is not surprising. One might, however, won- 
der why it works so much better for the ground state than 
for these excited states. A perhaps even more revealing 
question is why the data follows the GGRT energy spec- 
trum so closely even in the regime where the low energy 
expansion breaks down. 

The situation is similar, though slightly more intri- 
cate for excited states containing both a left- and a 
right-moving phonon each with one unit of longitudinal 
momentum. These two-particle states are conveniently 
characterized according to their transformation proper- 
ties under the unbroken 50(2). They break up into one 
scalar, one pseudoscalar (or cquivalently an antisymmet- 
ric tensor), and two components of a symmetric trace- 
less tensor. The effective string theory predictions are 
again universal up to £g/R 5 in the £ s /i?-expansion. How- 
ever, for these levels the universal terms of order £^/R 5 
for a relativistic string differ from those in the GGRT 
theory 7J. In the GGRT-theory, all these two-particle 
states are degenerate. For the relativistic string, repre- 




FIG. 2: This plot shows AE = E - Rji\ as a function of 
the length of the flux tube for the lowest lying states contain- 
ing both left- and right-movers. The scalar and pseudoscalar 
states are shown in blue and red, respectively. The spin-2 
states are shown in green. The solid lines show the theoretical 
predictions derived in the remainder of the paper, including 
the worldsheet axion in addition to the Nambu-Goto fields. 
The dotted lines show the prediction of the ^ s /i?-expansion. 
The gray dashed line is the GGRT prediction. 



sentations with different spins are split. The splitting 
originates from the Polchinski-Strominger (PS) interac- 
tion [7H9] and is proportional to D — 26. The numerical 
results of [5] are shown in Fig. [2] along with the pre- 
dictions of the GGRT theory, the derivative expansion 
of the relativistic theory, and the theory predictions we 
will discuss later. The scalar and symmetric tensor lev- 
els follow the GGRT prediction rather closely at large 
radius, but for small radii a splitting is clearly visible. 
The splitting qualitatively agrees with the one predicted 
by the universal PS terms. However, the ^ s /i?-expansion 
breaks down at radii that are so large that a quantitative 
comparison is impossible. Even more noticeable is that 
the pseudoscalar state strongly deviates from the GGRT 
model. 3 This data clearly calls for an alternative to the 
standard £ S /R expansion. 

Let us first identify the physical reason for the failure 
of the £ S /R expansion for excited states. To this end, it 
is instructive to first inspect the properties of the £ S /R 
expansion in the GGRT theory itself. Its exact spectrum 



3 One also notices a splitting between the two components of the 
symmetric tensor, but this splitting is due to lattice effects [2]. 



3 



E, 



GGRT 




Att 2 (N -Nf 

b: 2 



in 



N + N-- 
6 



Here N and N are the levels of an excited string state 
counting the amount of the longitudinal momentum car- 
ried by the left- and right-moving phonons separately, so 
that 2ir(N — N)/R is the total longitudinal momentum 
of the state. This formula immediately reveals a techni- 
cal reason for the breakdown of the £ S /R expansion. For 
excited states the expression under the square root not 
only involves powers of £ S /R, but also contains factors 
of 27T N, which turn the £ S /R expansion into a diverging 
asymptotic series even for relatively large values of R. 

To find a remedy let us reformulate the issue in more 
physical terms. In general, the energy of an excited state 
is of the form 

E = e- 1 £(p i £ a ,t s /R) 

where pi are the momenta of individual phonons com- 
prising the state. These are quantized in units of 2ir/R 
in a free theory, but we will see that this is modified in an 
interacting theory. So in reality the naive £ S /R expansion 
is a combination of two physically different expansions. 
The first expansion is an expansion in the softness of in- 
dividual quanta compared to the string scale, i.e. in pi£ s . 
The second expansion is a large volume expansion, i.e. 
an expansion in £ S /R. 

The key idea of our approach is to improve the conver- 
gence by disentangling the two expansions. To achieve 
this, we divide the calculation of the finite volume spec- 
trum into two steps. We first calculate the (infinite vol- 
ume) S-matrix. This can be done perturbatively pro- 
vided the center of mass energy of the colliding phonons 
is small in string units. We will call this the momentum 
expansion. We then calculate the finite volume energies 
from this S-matrix. Conceptually, this step is not re- 
lated to the perturbative expansion used to calculate the 
S-matrix. However, a prescription for this step for a 
given S-matrix only exists in two cases. First, there is a 
perturbative procedure due to Liischer [10l [11] which is 
routinely used in lattice calculations and applicable for 
states with energies below the inelastic threshold. Sec- 
ond, for two-dimensional integrable theories, there is the 
Thermodynamic Bethe Ansatz (TBA) [I2"1 IT5]. a set of 
integral equations, whose solution yields the exact spec- 
trum of the theory on a circle from the S'-matrix. 

In our case the particles are massless and we cannot 
rely on Liischer's method. However, let us ignore this for 
now and give the prescription to extract the finite volume 
spectrum from the S matrix in a way that allows us to 
discuss both methods in parallel. For massless particles, 
it is convenient to choose pi to be positive and to divide 
particles into left- and right-movers. The energy of the 



state then takes the form 

E = J2Pi+W E 
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where We represents winding effects from virtual quanta 
traveling around the circle. The particle momenta satisfy 
a modified periodicity condition, 



e(Pi) = 2-Krii/R 



where 
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and rii are positive integers. Interactions thus modify the 
quantization condition for momenta in two ways. First, 
real particles scatter with each other explaining the in- 
finite volume scattering phase shift 25(pi,pj) in Eq. (j3j. 
Second, it is modified by winding corrections represented 
by W P . 

In Liischer's approach the winding corrections W E , 
Wp are calculated perturbatively. For theories with a 
mass gap /x they are exponentially suppressed as e~^ R . 
In analyses of lattice results it is thus common to use the 
massive version of equations 0, ([2]), with We,p = 0. 
In the context of integrable field theories, the resulting 
simplified equations are known as the asymptotic Bethe 
Ansatz. 

In massless theories winding corrections are only 
power-law suppressed, requiring us to work harder and 
to use insights from integrable theories. For a generic 
S'-matrix, even writing the excited state TBA equations 
in closed form can be challenging, and a certain amount 
of guesswork is typically involved. However, different ex- 
amples (e.g., O rjS] ) strongly suggest that for purely 
elastic scattering the excited state TBA equations take a 
rather simple universal form 
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where the /(r) subscript refers to left(right)-movers, and 
the densities fu r )(p) are reminiscent of a thermal bath at 
temperature 1/i?, 
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The pseudo-energies ei( r )(p) are then determined by solv- 
ing the integral TBA equations 

In the GGRT theory the phase shift takes the simple 
form 
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allowing for an exact solution of the TBA equations [TS] . 

The worldsheet theory of a QCD flux tube is not inte- 
grable and the situation seems hopeless. However, as a 
consequence of the shift symmetry the theory is weakly 
coupled at low energies and the low energy scattering 
is dominated by purely elastic processes. Furthermore, 
the low-energy scattering amplitudes agree with those of 
the GGRT theory, which is intcgrable and for which we 
do know the exact form of TBA equations and winding 
corrections. We thus claim that we can use the TBA 
equations for the GGRT theory as a zeroth order ap- 
proximation and incorporate higher order contributions 
in the momentum expansion as corrections to the scat- 
tering phase shift. In principle, this should be done both 
in asymptotic and winding parts of the TBA equations. 
In practice, we include corrections only in the asymptotic 
part and use the GGRT phase shift for windings. This is 
justified by the form of Eq. It ensures that the wind- 
ing corrections receive their dominant contributions from 
virtual quanta with momenta below l/R, much softer 
than the real quanta. It is also important to note that 
the momenta of the real quanta for multiparticle states 
containing both left- and right-moving phonons are softer 
than the free theory estimate 2irrii/R. This is immedi- 
ate from pj),(|3]) and the fact that the GGRT phase shift 
^ grows in the UV. This indicates that one should ex- 
pect a better than naive agreement between the QCD 
flux tube and the GGRT spectrum. As soon as the rel- 
evant momenta are small enough, the two theories have 
similar infinite volume S'-matrices and as a result should 
have similar finite volume spectra. This fact is lost in the 
conventional perturbative expansion. 

Let us first apply this logic to the purely left-moving 
states. For these states the asymptotic Bethe Ansatz is 
trivial. The GGRT winding corrections are small, and 
one expects the spectrum to be close to that of a free 
theory. The dotted line in Fig. [T] shows that this expec- 
tation is correct. This eliminates the mystery for these 
states. 

For states containing both left- and right-movers, we 
need to take into account corrections to the GGRT phase. 
The leading one-loop correction to the amplitude is uni- 
versal and takes the PS form [7] 

llf 4 

2Sp,s(pi,p r ) = ±—^( PlPr ) 2 , (7) 
LZn 

where "+" refers to the scalar and pseudoscalar channels 
and " — " to the symmetric tensor channel. 

Including this contribution in the asymptotic Bethe 
Ansatz results in a significant improvement for the scalar 
and symmetric tensor levels. The improved theoretical 
control makes it manifest that the anomalous behavior 
of the pseudoscalar level cannot be blamed on the break- 
down of the perturbative expansion and a qualitatively 
new ingredient is needed. The energy of the anomalous 
level is roughly independent of the radius. This suggests 
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FIG. 3: This plot shows the scattering phase shift S for two 
Goldstone bosons as a function of the center of mass momen- 
tum in the symmetric traceless, scalar, and antisymmetric 
channel in the top, middle, and bottom panel, respectively. 
The solid and the long dashed lines show the theoretical pre- 
diction with and without the worldsheet axion, respectively. 



that the most straightforward way to explain this level 
is the introduction of a massive pseudoscalar particle <fi 
on the worldsheet. The leading interaction compatible 
with non-linearly realized Lorentz invariance for such a 
particle is a coupling to the topological invariant known 
as the self-intersection number of the string worldsheet 

S mt = £ J d^K^Kpe^e^ , (8) 

where K % ai is the extrinsic curvature of the worldsheet. 
The existence of this worldsheet #-term for a string 
in a four-dimensional target space was pointed out by 
Polyakov [TB], and it was suggested that it should be 
generated on the flux tube worldsheet in the presence of 
the bulk 6*-term [17] . Given this coupling, it is natural to 
refer to the field <p as the worldsheet axion. 

The axion appears as a resonance in the scattering of 
Goldstone bosons with antisymmetric flavor wave func- 
tion and it also contributes to the scattering in the scalar 
and symmetric tensor channels through t- and it-channel 
diagrams. It is thus readily included in the TBA equa- 
tions. By fitting the two free parameters (the axion mass 
m and the coupling a) to the data, we obtained the spec- 
trum shown as solid lines 4 in Fig. [2] which corresponds to 



4 The lines are dashed where |7]l becomes comparable to |6j} and 
our approximations are unreliable. 
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m£ s — 1.85^q'q3, confirming the heuristic analysis of [2J, 
and a = 9.6 ± 0.1. The error bars represent the statisti- 
cal uncertainty only. Based on a comparison of the two 
symmetric tensor levels and a comparison of the states 
with zero and one unit of longitudinal momentum, we 
estimate the systematic and theoretical uncertainties to 
be about a factor of five larger. 

Further support for existence of this axion comes from 
data for the next excited level in the pseudoscalar chan- 
nel. We reverse the logic and use the TBA equations to 
determine the scattering phase shift from the finite vol- 
ume spectrum, which is the standard approach in lattice 
QCD. The resulting phase shifts for the pseudoscalar, 
scalar and symmetric tensor channel are shown in Fig. [3] 
For the pseudoscalar, it exhibits a characteristic reso- 
nances shape. The phase shift extracted from the data 
for the pseudoscalar state we discussed so far is shown in 
light red. The dark red points show the phase shift ex- 
tracted from the data for the next excited pseudoscalar 
state also taken from [5J. 

In summary, the TBA approach provides better theo- 
retical control over flux tube spectra than the standard 
perturbative expansion. We conclude that existing lat- 
tice data provides strong evidence for the existence of a 
new particle - the worldsheet axion. In a forthcoming 
publication [TS] we will present the details of our analy- 
sis and elaborate on the diagrammatic interpretation of 
the TBA method. 
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